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MAXIMUM PRINCIPLES 
FOR THE RELATIVISTIC HEAT EQUATION 

EVAN MILLER AND ARI STERN 


Abstract. The classical heat equation is incompatible with relativity, 
since the strong maximum principle allows for disturbances to propagate 
instantaneously. Some authors have proposed limiting the propagation 
speed by adding a linear hyperbolic correction term, but then even a weak 
maximum principle fails to hold. We stndy a more recently introduced 
relativistic heat equation, which replaces the Laplace operator by a 
quasilinear elliptic operator, and show that strong and weak maximum 
principles hold for stationary and time-varying solutions, respectively, 
as well as for sub- and supersolutions. Moreover, by transforming the 
equation into an equivalent form, related to the mean curvature operator, 
we prove even stronger tangency and comparison principles. 


1. Introduction 

1.1. Background. It is well known that the classical heat equation, 

(1) ut = An, 

allows for disturbances to propagate with infinite speed, in violation of special 
relativity. This is a direct consequence of the strong maximum principle, 
which states that if u attains an interior maximum at some positive time, 
then it must be constant for all previous times. This implies that, if a 
constant solution is perturbed locally, then the perturbation is instantly 
detected at arbitrarily distant points. 

To model heat conduction relativistically, various authors have suggested 
replacing the parabolic heat equation by a linear hyperbolic equation, such 
as the telegraph eguation, 

(2) c~‘^utt + ut = Au, 

where the constant c > 0 denotes the speed of light. Indeed, this equation 
obeys the relativistic “speed limit” and reduces to 0 in the limit as c —)• oo. 
(For example, see Gurtin and Pipkin [^.) However, as a hyperbolic equation, 
it differs from 0 in at least two important respects. First, it lacks the 
regularity and smoothing properties of 0 : solutions of 0 need not be 
smooth and can even contain “thermal shock waves” (cf. Tzou [I^)- Second, 
it does not even satisfy a weak maximum principle: for example, 0 allows 
for two “heat waves” traveling towards one another to combine into a larger 
wave, in violation of this principle. 
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Brenier [2] introduced an alternative approach to relativistic heat conduc¬ 
tion, based on optimal transportation theory. This results in the quasilinear 
parabolic equation 

uDu 

ut = div — , . 

+ c~‘^\Du\^ 

which we call the relativistic heat equation. If > 0, we again obtain 0 in 
the limit as c —)• oo. (By contrast, if rt < 0, the limit is instead the ill-posed 
backward heat equation, ut = —Au.) Brenier’s starting point was to observe 
that Q can be understood as gradient descent of the Boltzmann entropy 
functional f u log u dx, where the gradient is with respect to the Wasserstein 
metric, corresponding to optimal transportation with the nonrelativistic 
kinetic energy cost function k{v) = \\v\‘^ (cf. Jordan et al. [7j). Using instead 
the relativistic cost function 


k{v) 



f I < c. 


oo, 


f I > c. 


yields ([^ as the gradient descent equation for Boltzmann entropy. For the 
remainder of this paper, we set c = 1. 

It is natural to ask whether the relativistic heat equation 0 satisfies a 
weak maximum principle, similar to that satisfied by 0 but not by ([^. The 
purpose of the present paper is to answer this question in the affirmative, 
and to give some related results on maximum principles for the relativistic 
heat equation. 


1.2. Outline of the paper. First, in Section 2 we consider stationary 


solutions to the relativistic heat equation, which we call relativistically har¬ 
monic functions (by analogy with harmonic functions, which are stationary 
solutions to the classical heat equation); we also consider subsolutions and 
supersolutions, which we call relativistically subharmonic and superharmonic, 
respectively. A crucial component of our analysis is Lemma 2.2, which 


transforms the quasilinear elliptic operator in Q into a more convenient 
form, related to the mean curvature operator. Using this transformed formu¬ 
lation, we prove a strong maximum (minimum) principle for subsolutions 
(supersolutions), as well as even stronger tangency and comparison principles. 
Finally, we prove that relativistically harmonic functions are real analytic, 
and use this to give another, more elementary proof of the strong maximum 
principle. 

Next, in Section 3, we consider time-dependent solutions of the relativistic 
heat equation, along with subsolutions and supersolutions. While finite prop¬ 
agation speed (i.e., relativity) precludes the possibility of a strong maximum 
or minimum principle, much less an even stronger tangency principle, we 
show that comparison and weak maximum/minumum principles do hold. 
Finally, we discuss one possible direction for future work, in which a stronger 
maximum/minimum principle and tangency principle might be shown to 
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hold on light cones, which would still be consistent with hnite propagation 
speed. 


2. The elliptic case 

2.1. Solutions, subsolutions, and snpersolutions. 


Definition 2.1. Given an open set U C 
tt > 0 is relativistically subharmonic if 

u Du 


( 4 ) 


div 


y/u^~+~\Du\^ 


\ a function u G C'^{U) with 


> 0 , 


relativistically superharmonic if 
(5) div 


u Du 


V^u2 + |Du|2 


and relativistically harmonic if 
(6) div 


u Du 


< 0 , 


= 0 , 


y/u^ + \Du\^ 
i.e., if it satisfies both Q and ([^. 

For the purposes of the subsequent analysis, these expressions are less 
than ideal. In addition to the u > 0 restriction, the coefficients depend on 
both u and Du, while several results on quasilinear elliptic operators require 
dependence on Du alone. However, by making the substitution w = log u, 
we now show that it is possible to obtain an equivalent formulation that is 
valid for all real-valued w, and where the coefficients depend only on Dw, 
not on w itself. 

Lemma 2.2. Given u G C‘^{U) with u > 0, let w = logtt, and define the 
quasilinear operator 

D‘^w{Dw, Dw) 


Qw = Aw — 


+ \Dw\^ 


l + \Dw\^ 

Then u is relativistically subharmonic if and only if 

(7) Qw > 0, 

superharmonic if and only if 

(8) Qw < 0, 

and harmonic if and only if 

(9) Qw = 0. 

Proof. Since u > 0, observe that 


u Du 


u 


Du 


e^Dw 


y/u'^ p |Du|2 y/YfifiDw^ 

Y ' 
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Taking the divergence of this expression, a short calculation gives 


div 


e^Dw 


[Aw + l-Dtcp) (l + \Dw\‘^) — D‘^w{Dw, Dw) 


^l + \Dw\^ 


(1 + 

Substituting this into each of Q-Q, dividing by e"'(l + \Dw\‘^) 
and rearranging yields respectively. 

Lemma 2.3. The quasilinear operator Q is non-uniformly elliptie. 


> 0 , 

□ 


Remark 2.4. This essentially amounts to the non-uniform ellipticity of the 
well-studied mean curvature operator ‘Aflw = [l + \Dw\‘^^ Aw — D‘^w[Dw, Dw) 
(cf. Gilbarg and Trudinger [3]), but the demonstration is sufficiently brief 
that we give it here. 


Proof. The principal part of Q can be written as {Dw)DiDjW, where 


ab(p) = 


pipj 

1 + IpP ’ 


with <5*-^ denoting the Kronecker delta. The matrix p^pp is symmetric with 
rank 1, so its only nonzero eigenvalue (having multiplicity 1) is its trace, \p\^. 
Diagonalizing, it follows that a^^{p) has eigenvalues 

^ l + |p|2’ A2 = --- = A„ = 1, 

which are positive; however, Ai is not bounded uniformly away from zero. □ 


2.2. Maximum/minimum, tangency, and comparison principles. We 

begin by giving the strong maximum and minimum principles for subsolutions 
and supersolutions, respectively. 


Theorem 2.5. Suppose U C M"" is open and eonnected. 

(i) If w £ C‘^{U) satisfies Qw > 0 and attains an interior maximum in 
U, then w is constant. 

(a) If w £ C‘^[U) satisfies Qw < 0 and attains an interior minimum in 
U, then w is constant. 

Consequently, if u £ C‘^[U) is relativistically subharmonic (superharmonic), 
then it satisfies the corresponding strong maximum (minimum) principle. 


Proof. Using Lemma 2.3 and the fact that the elliptic operator contains no 
zeroth-order terms in w, the statements for w follow immediately by applying 
Hopf’s strong maximum principle. (On the applicability of Hopf’s principle 
to nonlinear elliptic inequalities by “freezing” coefficients, which is lesser- 
known than the linear case, see Pucci and Serrin (9j.) The corresponding 
statement for u then follows by the monotonicity of w e^ = u. □ 


In fact, using the fact that Q is independent of tc, we can establish 
an even stronger “tangency” principle, which implies the preceding strong 
maximum/minimum principle as a special case. 
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Theorem 2.6. Suppose w,w' G C'^{U) satisfy Qw > Qw' and w < w' in 
U. If w = w' at some point x ^ U, then w = w' in U. Consequently, if 
u, u' G C'^{U) satisfy 


( 10 ) 


div 


u Du 


> div 


u' Du' 


^u^ + \Du\^ ~ ^{u'Y + \Du'\^^ 


and u < u' in U, and if u = u' at some point x £ U, then u = u' in U. 


Proof. Since Q is elliptic, and since its coefficients are independent of w and 
continuously differentiable (in fact, analytic) in Dw, the result is obtained 
by applying the tangency principle for nonlinear elliptic operators in Pucci 
and Serrin u Theorem 2.1.3]. Finally, as in the proof of [Theorem 2.5l the 


corresponding statement for u and u' follows by the monotonicity of the 
exponential function. □ 


Remark 2.7. Note that Theorem 2.5 is a special case of Theorem 2.6, where 
we take either w' = M, 
minimum attained by w'. 


the maximum attained hy w, or w = m', the 


Finally, we establish a comparison principle, which relates the values of w 
and w' on the boundary with those in the interior. Taking either w or w' to 
be constant, we get the weak maximum principle as a special case. 


Theorem 2.8. Suppose w,w' G C‘^{U) n CiU) satisfy Qw > Qw' in U. If 
w < w' on dU, then w < w' in U. Consequently, if u,u' G C‘^{U) n C{U) 
satisfy the inequality (10) inU and u <u' on dU, then u <u' in U. 


Proof. Again, since Q is elliptic, and since its coefficients are independent of 
w and continuously differentiable (in fact, analytic) in Dw, we may apply 
the nonlinear elliptic comparison principle in Pucci and Serrin 0 Theorem 
2.1.4] or its quasilinear counterpart in Gilbarg and Trudinger [H, Theorem 
10.1]. As above, the corresponding statement for u and u' follows by the 
monotonicity of the exponential function. □ 


2.3. Analyticity and an elementary proof of the strong maximum 
principle for relativistically harmonic functions. Finally, we show 
that relativistically harmonic functions are analytic, and we use this to give 
an elementary, self-contained proof of the strong maximum principle for 
solutions (but not sub- or supersolutions) using analyticity rather than the 
machinery of the Hopf principle. 

Theorem 2.9. Every solution w G C‘^{U) of Qw = 0 is real analytic. 
Consequently, every relativistically harmonic function u is real analytic. 


Proof. Analyticity of w follows from [Lemma 2.3 by classical elliptic theory 
(e.g., Hopf [U|, Morrey 0), since the coefficients depend analytically on Dw. 
Analyticity of u then follows immediately from Lemma 2.2 □ 


We now give an alternative proof that, if a relativistically harmonic 
function attains an interior maximum, then it must be constant. (The proof 














6 


EVAN MILLER AND ARI STERN 


of the minimum principle is essentially identical, modulo the direction of the 
corresponding inequalities.) 


Proof. Suppose u attains an interior maximum u{x) = M at some x £ U. 
Since u is real analytic by Theorem 2.9, the maximum at x must be isolated 
unless u is constant. In the latter case, we’re done. Otherwise, the maximum 
is isolated, so there exists an open ball B{x, r) C U oi radius r > 0 such that 
Du{y) ■ {y — x) <0 for all y G B{x,r). That is, Du • i/ < 0 on dB{x, s) for all 
0 < s < r, where r'(y) = ^(y — x) is the outer unit normal at y G dB{x, s). 
Now, since 

u Du 

= 0 , 


div 


^yu?~+~\Du^^ 


integrating over B(x,s) and applying the divergence theorem implies 


f " 

JdB(x,s) \/u^ P \Du\^ 


Du ■ V d5 = 0. 


However, since —p=d^=- > 0 and Du • z/ < 0 on dB(x, s), it follows that 
^fu^ + \Du\^ 

Du ■ v = 0 on dB{x, s). Since this holds for all 0 < s < r, the function u is 
constant along radii of B(x,r), and hence constant on B{x,r). Explicitly, if 
y G B{x,r) and u = v{y), then 


u{y) 



d 

ds 


u{x + vs) ds 



Du{x + vs) • ds = 0, 


so u = M on B{x,r). Hence, the nonempty and relatively closed set 
u~^[{M}) nU is also open, so the result follows by the assumption that U 
is connected. □ 


3. The parabolic case 


3.1. Solutions, subsolutions, and supersolntions. We now turn our 
attention to time-dependent solutions of the relativistic heat equation, along 
with subsolutions and supersolutions. Throughout this section, we denote 
the spacetime domain by {x,t) G Ut x (0,T], where U C M” is an open set 
and T > 0, and the parabolic boundary of Ut by T-p = Ut — Ut- We say 
that u G C\{Ut) if u{x,t) is in x and in t for all {x,t) £ Ut- (This is 
consistent with the notation found, e.g., in Evans [3].) 


Definition 3.1. Given u G Ci{Ut) with u > 0, we say that n is a subsolution 
of the relativistic heat equation if 


a supersolution if 


Ut — div 


uDu 

■\/u^ + \Du\^ 


Ut — div 


u Du 

■\/u‘^ + \Du\‘^ 


< 0 , 

> 0 , 
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and a solution if 


ut — div 


uDu 


+ \Du\^ 


= 0 , 


i.e., if u is both a subsolution and supersolution. 


As in the previous section, we use the change of variables w = log u to 
transform the problems above into a more convenient form, where the elliptic 
coefficients depend only on Dw rather than on w itself. 

Lemma 3.2. Given u G Ci{Ut) with u > 0, let w = logu, and define the 
quasilinear operator 

^l + \Dw\^ 

Then u is a subsolution of the relativistic heat equation if and only if 

wt — Qw < 0, 

a supersolution if and only if 

Wt — Qw > 0, 

and a solution if and only if 

Wt — Qw = 0. 


Proof. Since w = log u, we have 


Ut jjj 

Ut = u— = e wt, 
u 


so the result follows by the calculation given in the proof of Lemma 2.2 □ 


3.2. Comparison and weak maximum/minimum principles. The rel¬ 
ativistic heat equation, as its name suggests, does not permit instantaneous 
propagation of disturbances (cf. Andreu et al. [T]), and hence we cannot 
expect a strong maximum principle or tangency principle to hold in the 
time-dependent case. However, we now show that a comparison principle 
does still hold, which implies a weak maximum/minimum principle as an 
immediate corollary. 

Theorem 3.3. Suppose w,w' G CKUt) Pi C{Ut) satisfy 

Wt — Qw < w't — Qw' 

in Ut- If w < w' on Ft, then w < w' in Ut- Consequently, if u,u' G 
CI{Ut) n C{Ut) satisfy 

u Du , u' Du' 

Ut — div — , = < Ut — div — , 

^/uf^jDu^ ^yiu')2 + \Du'\^ 

in Ut, and if u < u' on Ft, then u < u' in Ut- 
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Proof. The proof is essentially a parabolic adaptation of the quasilinear 
elliptic comparison principle (Gilbarg and Trudinger [U Theorem 10.1]). 
First, we rearrange the inequality to obtain 

{w — w')t — {Qw — Qw') < 0. 


Next, 

Qw — Qw' = {Dw)DiDj{w — w') 

+ {Dw) — {Dw')\DiD jw' + b{Dw) — b{Dw'). 

Letting z = w — w', we then define the linear operator Lz = t)DiDjZ + 
P{x, t)DiZ so that 

{x, t) = '(P^{Dw), 

b\x,t)DiZ = [P^Dw)-P^{Dw')]DiDjw' + biDw) -b{Dw'). 

(In the last step, it is crucial that a*-^ and 5* are independent of w, w' .) Hence, 
Qw — Qw' = Lz, so we have 


zt — Lz < 0 


in Ut, and 2 < 0 on Ty. By the parabolic weak maximum principle (cf. Evans 
[21 Section 7.1, Theorem 8]), we conclude that z < 0 in Lr, i.e., w < w', 
which completes the proof. 

Finally, as in the previous section, the corresponding statement for u and 
u' follows by the monotonicity of the exponential function. □ 

Corollary 3.4. Let we Cl{UT)r\C{UT). 

(i) If wt — Qw < 0, then maxjj^ w = maxp^ w. 

(ii) If Wt — Qw > 0, then min^^ w = minrj, w. 

The corresponding statements hold for subsolutions and supersolutions of the 
relativistic heat equation. 


We also mention another corollary of Theorem 3.3, which establishes 


monotonicity and uniqueness properties for solutions to the relativistic heat 
equation 


Corollary 3.5. Suppose w,w' e Cf{UT)nC{UT) are both solutions to the 
relativistic heat equation. Ifw<w' on F^, then w < w' in Ut- Furthermore, 
if w = w' on Tt, then w = w' in Ut. 


Proof. Since w and w' are solutions, wt — Qw = w't — Qw' = 0. Hence, the first 
statement follows immediately by applying [Theorem 3.3 , while the second 
statement follows by applying it again with w and w' interchanged. □ 
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3.3. Conjectured strong maximum/minimum and tangency prin¬ 
ciples in light cones. We conclude with a brief discussion of possible 
directions in which the results of this section might be strengthened in the 
future—but which would also require the introduction of new techniques. 

The strong maximum (minimum) principle for the classical heat equation 
says that, if a subsolution (supersolution) u attains an interior maximum 
(minimum) at a point {x, t) G Ut, then u must be constant on Ut = U x (0, t]. 
As mentioned above, we cannot hope for this to hold for the relativistic heat 
equation—at least not for the cylinder Ut- However, it seems likely that an 
analogous statement might hold on the backwards light cone 

U(x,t) = {(^,'r) ^Ut :\x-^\<t-T}, 

which contains only points that can affect (x, t) without violating relativity. 
(For c —)• oo, this approaches the cylinder Ut-) The following conjecture 
states a version of this principle restricted to this backwards light cone. 

Conjecture 3.6. If w G Ci{Ut) H C{Ut) satisfies wt — Qw < 0 and attains 
an interior maximum M at a point (x, t) G Ut, then w = M in the backwards 
light cone U(^x,t)- Likewise, if w satisfies wt — Qw > 0 and attains an interior 
minimum m at a point (x,t) G Ut, then w = m in 

We also expect that an even stronger tangency principle might hold on 
the backwards light cone, as follows. 

Conjecture 3.7. Suppose w,w' G C‘1{Ut) H C{Ut) satisfy 

Wt — Qw < w't — Qw' 

and w < w' in the backwards light cone U(^x,t) of a point (x,t) G Ut- If 
w{x,t) = w'{x,t), then w = w' in 
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